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_ METHOD AND APPARATUS FOR SCHEDULING SWITCHED 

;g MULT I BEAM ANTENNAS IN A MULTIPLE ACCESS ENVIRONMENT 

si i 

ill Related Applications 

III This application is based on provisional patent 

].5 application No. 60/155, 736 filed September 23, 1999 entitled 

,g "Scheduling of Switched Multibeam Antennas in a Multiple 

m Access Environment", the disclosure of which is incorporated 
herein by reference. 

Field of the Invention 

10 The invention pertains to intelligent switched-beam smart 

antennas for radio communications systems with multiple 
transmitters per receiver. The invention is particularly 
suitable for cellular communications, but is applicable to any 
other type of wireless communication, such as radar, sonar, 

15 wireless local area networks, wireless local loops, and 
satellite communications systems. 
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Background of the Invention 

Beamf orming techniques are well known in radio 
communications systems in which a single receiving station 
must be capable of receiving signals from multiple signal 
5 sources, i.e., transmitters (or a single source, if that 
source is moving) for increasing the focus or gain of the 
receiver in the direction of the signal (s). In particular, a 
plurality of receiving elements, i.e., antennas, can be 
focused to receive signals in a particular angular portion of 
ig6 the receiving area while rejecting other signals and 
Jl: interference outside of that angular portion of the receiving 

area. This can be done by using multiple, directional 
!*; antennas focused on different angular portions of the 

geographic area of reception surrounding the antenna array. 
jK? Alternately, the electrical outputs of an array of omni- 
Ji: directional receiving antennas can be electronically weighted 

w and combined in various manners so as to create a plurality of 

beam signals, each beam being particularly focused on a 
particular angular portion of the geographic region 
20 surrounding the antenna array. The output signals of the 

antennas, of course, are radio frequency (RF) signals. These 
outputs may or may not be frequency downconverted to an 
intermediate frequency (IF) prior to any processing. Thus, 
the circuitry for combining the antenna output signals 
25 operates in the RF or IF band. 



PATENT 3 Docket No. 2 4,338-A USA 

Beamf orming is commonly used in wireless communication 
systems, such as in cellular telephone system base stations, 
for increasing reception performance. For instance, by 
focusing a receiving beam on a particular transmitter, that 
5 transmitter need transmit with much less power to be 

adequately received by the base station relative to a system 
using a base station with only a single omnidirectional 
antenna and no beamf orming. 

There are three general beamforming techniques in common 
i=i> use at this time, namely, fixed beamforming, adaptive 
j£f beamforming and switched beamforming. Fixed beamforming can 
; 2 be implemented with an array of multiple, directional antennas 

or by use of an array of multiple, omnidirectional antennas 
!L. with electronic beamforming circuitry. In fixed beam systems, 
5t$ the beams are permanent and fixed. 

In adaptive beamforming systems, the beamforming is 
performed electronically, rather than physically with 
directional antennas. Further, the antenna electrical RF 
output signals are combined using adaptive combining circuitry 
20 such that a receiving beam can track a moving signal source. 

Finally, in switched beam systems, the RF electrical 
output signals of multiple, omnidirectional antennas of an 
antenna array are combined by RF combining circuitry to 
generate a plurality of fixed beam output signals, each 
25 covering an angular portion of the geographic area surrounding 
the array. However, in switched beam beamforming techniques, 
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an algorithm is employed to select the one of the fixed beam 
output signals that provides the best reception of the signal 
source (s) at any given instant in the region (the region being 
termed a cell in cellular communications systems) . As 
5 conditions change, the selected beam also changes accordingly. 
The switching rates typically used in switched beam 
beamforming systems are typically substantially slower than 
the data rate of the system. One common scheme is to switch 
about once a frame or even slower, a frame typically 

jess. 

3r€) comprising several hundreds of data symbols. 

i Hi 



Switched beamforming systems are commonly in use today 
.jf with base stations of cellular telephone communications 

■*. s systems, including code division multiple access 

iws: 
s 

;L, (CDMA) cellular telephone communications systems. 

JS Figure 1 is a block diagram of a switched beam 

IZ beamforming system in accordance with the prior art. An 
•~ antenna array 12 comprising multiple antennas 12 x , 12 2 , 12 N 

are arranged at a receiving station 10, e.g., at the top of a 
cellular telephone system cell tower. The outputs of those 
20 antennas are fed into an analog NxB beamforming circuit 14, 
wherein N is the number of antennas and B is the number of 
beams generated by the beamforming circuitry from the antenna 
signals. The analog beamforming circuitry 14 weighs and 
combines the RF antenna output signals on lines 13 x , 13 2 ,..., 
25 13 N in accordance with a scheme dictated by beamforming 

control circuit 21 to produce the beam signals 15 lf 15 2 ,..., 
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15 B , each beam focused on an angular portion of the reception 
area. The number of beams, B, typically is less than or equal 
to the number of antennas, N, in the antenna array. Each beam 
signal 15 a/ 15 2 ,..., 15 B is passed through frequency down 
converting circuitry 16 T , 16 2 ,..., 16 B for converting the beam 
signals from the RF or IF frequency range to the baseband 
range. Each of those signals is digitized by an analog to 
digital (A/D) converter 18 x , 18 2 ,..., 18 B . The outputs of the 
A/D converters are then each input to a fading multipath and 
multi-user channel estimation circuit, 20!, 20 2 ,..., 20 B . Each 
of those circuits generates L output signals for each of M 
simultaneous transmitters at use in the given geographic area, 
where L is the number of paths per transmitter that the 
circuitry is designed to process simultaneously (e.g., 
typically around 3 or 4) and M is the number of simultaneous 
transmitters using the receive station. Those outputs are 
input to a minimum variance selector 22 for all of the paths 
and users. The minimum variance selector generates a path 
estimate and a path estimate error for each path of each user, 
i.e., L x M path estimates and L x M path estimate errors. 

It should be understood by those skilled in the related 
arts that all circuitry subsequent to the analog-to-digital 



converters, i.e., circuits 20^ 20 2 ,..., 20 B and 22 typically 
would digital and may be implemented by a digital signal 
processor (DSP) , a processor, a microprocessor an application 
specific integrated circuit, a finite state machine, a 
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programmed general purpo/se computer or any other equivalent 
(hereinafter collectively "DSP" ) . In fact, because of the 
large amount of processing needed to generate the path 
estimates and path estimate errors from the beams, the flat 
fading multipath and multiuser estimation and the minimum 
variance selection typically is performed by a bank of DSPs, 

In the analog N x B beamforming circuitry 14, for each 
beam, dedicated RF or IF combining circuitry is necessary, 
i.e., there are N copies of essentially identical circuitry 
H) for processing the N antenna signals on lines 13 x , 13 2 ,..., 

13 B . RF or IF band hardware is much more expensive than DSPs 
^ and other baseband circuitry. Also, as shown in Figure 1, 
each beam requires a dedicated frequency down converting 
circuit 16 and a dedicated analog-to-digital converter 18. 
3f5 Accordingly, it is an object of the present invention to 

:~ provide an improved switched beam beamforming method and 

w apparatus . 

It is another object of the present invention to provide 
a switched beam beamforming apparatus which reduces the amount 
20 of RF or IF band circuitry needed. 

It is a further object of the present invention to 
provide a lower cost switched beam beamforming method and 
apparatus . 

It is yet another object of the present invention to 
25 provide a switched beamforming method and apparatus in which 
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the bulk of the signal processing requirements are transferred 
from the RF or IF band to the baseband. 

It is yet a further object of the present invention to 
provide a switched beam beamforming method and apparatus 
5 having equivalent or superior performance to prior art 

techniques utilizing much lower cost circuitry components. 

Summary of the Invention 

The invention is a switched beam beamforming method and 
ijf) apparatus for wireless communication receiving stations 
Ji: utilizing an array of antenna elements in which only a single 

beam is formed at any given instant and the beam is switched 
p at a very high rate, e.g., the chip rate of a direct sequence- 

!L code division multiple access communication system. An 
}J5 algorithm for scheduling the beamforming sequence is disclosed 
that optimizes performance by optimizing the 
signal/interference-plus-noise ratio for a given set of 
conditions at any instant in time. Those conditions include, 
for instance, the number and location of the various 
20 transmitters and the environmental conditions that dictate the 
various paths of the transmit signals in the reception area 
(e.g., cell). In particular, spatial diversity offered by 
antenna arrays for direct sequence-code division multiple 
access communication systems is exploited by an intelligent 
25 switched beam antenna at radio frequency level. The design is 
optimized to yield conditional mean estimates of the 
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communication channel during uplink transmission and compute 
minimum variance estimates of the communication channel by 
optimally combining the signals of the spatially distributed 
antennas at chip rate. 

Brief Description of the Drawings 

Figure 1 is a block diagram of an exemplary switched beam 
forming device of the prior art. 

Figure 2 is a block diagram of a switched beamforming 
device in accordance with the present invention. 

Detailed Description of the Invention 
I . INTRODUCTION 

One of the primary objects of the present invention is to 
provide a receiving circuit in which the received RF signals 
are immediately down converted to the base band and sampled 
after the antenna and the beamforming is performed in the 
baseband. 

In this specification, we design a new switched beam 
antenna array receiver that is particularly suited to direct 
sequence-code division multiple access (DS-CDMA) communication 
systems operating in multipath Rayleigh flat-fading channels. 
However, the invention is not limited to DS-CDMA and can be 
more broadly applied to other spread spectrum and multiple 
access communication systems, including, but not limited to, 
other CDMA communication methods and TDMA communication 
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methods. Also, the invention is not limited to Rayleigh flat- 
fading channels, but can be used in other channels, such as 
frequency selective channels . 

In accordance with the invention, the beamf orming 
technique selects the best beam from a set of fixed beams at a 
high frequency, such as chip rate, and generates a beam signal 
only for that selected beam at any given instant. Chip rate 
is the transition rate of the spreading sequence of the DS- 
CDMA scheme and is many times faster than the data rate of the 
data signal which is modulated with the spreading sequence in 
order to create the spread spectrum transmit signals. In this 
manner, when there are multiple simultaneous users (i.e., 
transmitters within the cell) , multiple beams can be selected 
during a single data symbol interval, thus providing maximized 
performance, e.g., the highest signal/noise-plus-interference 
ratio. A Kalman filter is used to estimate the attenuated 
spread spectrum signals for all users and all transmission 
paths . 

The theory, algorithms and operations that allow us to 
form only one beam at each instant in time with beam switching 
at a high rate, such as chip rate, will now be described in 
detail. The remainder of this specification is organized as 
set forth below. In section II, we describe the communication 
system, give the model assumptions and introduce our DS-CDMA 
switched-beam communication system as a linear, Gaussian, 
time-varying state-space model. In section II. 1, we formulate 
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the estimation objective, by seeking an optimal chip-rate 
beamf orming scheduling (switching) sequence that yields 
minimum variance estimates of the received spread spectrum 
signals. Section II. 2 presents some well known beamforming 
techniques. In section III, we present a method for computing 
the optimal scheduling for switched beam antenna arrays. 
Finally, we summarize our work in section IV. 
II . MODEL ASSUMPTIONS AND PROBLEM FORMULATION 

Consider an asynchronous DS-CDMA system with M transmitting 
users. The continuous-time baseband complex transmitted signal 
S m (t) of the mth user is given by 

S m W=^b m {f)c m (f)e J + (i) 

where P m , b m (t), c m (t) and t m denote the transmitted power, data 
symbol, the spreading signal, and phase of the mth user. The 
spreading signal c m (t) is given by 

00 

c m(t) = X V^^t " kT c> < 2 > 

k= -00 



where rect Tc (t) is a unit rectangular pulse of duration T c , and 
c k , m e {-1,1} denotes the kth chip of the mth user. The 
spreading code sequences {c k/m } for all users are known to the 
receiver . 

Let T b denote the symbol (or bit) duration. The ratio 
G=T b /T c is known as the processing gain. In this specification, 
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we assume G to be an integer for convenience. The data signal 
b m (t) is given by 

oo oo 

. b m (t)= X K„rect Tb (t - nT„)= £ b c km rect Tc {t - nT c ) (3) 

k=-oo k=-oo 



where b n , m e {-1,1} is the nth data bit of the mth user. Let 

5 b n,m = b 0+nG,m = b 1+ nG,m =• • • = b G-l+nG,m for n ^ { -~, ...,<*»} , m 6? 

!□ 

=0 {1,...,M}. Thus, we have subdivided the symbol b n/m into G 

I;] consecutive and equal symbols {b£ +nGfM , bJ +nGrm/ . . , b S-i +nG , m } • 

The period of eachb^ m is equal to the chip period T c . In this 
III specification, the data signals are treated as deterministic, 

! |J 

jg) but unknown, functions of time. 

Due to reflections from large objects, diffractions and 
scattering, multiple, delayed and attenuated copies of the 
same transmitted signals are induced on the different elements 
of the array. Here, in a preferred embodiment, we assume an N 

15 element uniform linear array. In an exemplary embodiment, the 
antenna elements are equally spaced along the x axis and the 
distance between two adjacent elements is d = A/2, where A 
denotes the wavelength of the carrier. 

The N vector channel impulse response from the mth user 

20 to the antenna array is given by 
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/=1 

where L m (t) is the known - possibly time varying - number of 
resolvable paths. 5(t) is the Dirac delta function. + m ,i(t) is 
the known direction of arrival (DOA) or angle of arrival (AOA) 
of the /nth user's signal via the Ith path. The AOA is 
measured anti-clockwise from the positive real axis. P m ,i(t), 
\j; m/1 (t) and T k/1 (t) are the signal attenuation, phase shift and 
time delay of the 2th multipath component for user m, 
respectively. a(cp) is known as the steering vector, and 
denotes the response of the antenna array to a signal 
impinging onto the array from an angle cp. For a uniform 
linear array, the steering vector is given by 

<*(♦) = [gi(<hai(<h / g2<<h a 2«h / / 9u(*) a N(*)] < 5 > 

where g n (+) denotes the directivity (gain) of the nth antenna 
element impinging on the array from angle 4». Here we assume 
for illustrative purposes that all the elements have unit 
directivity, i.e., g x (♦) = .. .=g N ( ♦) =1 . For our uniform linear 
array, we have 

a n «t>) = e Wn - 1)O0Btt) # Vn e {i n} (6) 

The receiving signal at the array elements is given by 
the superposition of all of the channel distorted users' 
signals, that is 
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M 

KO = Z s„, (0*^0,0 

= Z Z a i(f> mJ )a m l it)b m (/ - r m l (ty)c m (t - r mJ (/)) (7) 



where * denotes the convolution operator, and the complex 
amplitude oc m/1 (t) 

| <* m A0= JlP^p^iOe^^'^ (8) 

denotes the attenuated, phase shifted multiplicative channel 
fn interference at time t via the 2th path for the mth user, 

ii Prior to sampling, the signals from each element of the 

\Q antenna array are multiplied with a complex number and then 

|1| summed to form the array output. The process of combining the 

03) signals at the array is known as beamforming. The output of 
the beamforming array is given by 

y(t) = f H (t) r (t)+f H (t)w(t) (9) 
where w(t) is the vector representing the additive background 
noise, (*) H denotes the complex transpose operator, and f(t) 
15 denotes the time-varying beamforming vector. 

A discrete time version of Eq. (9) is given by sampling 
the output of the beamformer at fixed intervals. In this 
exemplary embodiment, the sampling rate is set to the chip 
rate. The scalar complex discrete-time beamformer output is 
20 thus given by 
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y ( kT c + A ) = f H ( kT c + A ) r ( kT c + A ) + f H ( kT c + A ) w ( kT c + A ) (10) 
for k=...,-2, -1, 0,1,2,... . The sampling offset is denoted 
by A, such that 0^A<T c . 
Model and Statistical Assumptions: 

We now present some important assumptions for our 
analysis. Those of skill in the art will recognize that other 
assumptions are possible while remaining within the scope of 
the present invention. 

Flat fading : We model the communication channel as flat- 
fading. The delay spread, which is defined as the inverse of 
the channel coherence bandwidth, is assumed much smaller than 
the symbol period. Here, we further assume that, 0<x m/1 ( t ) <T C , 
for all times t, users m, and propagation paths 1. Let the 
maximum number of paths L m (t) for each user to be equal to L. 
Sampling : We chose our sampling instances, such that the 
sampling offset A satisfies the following condition 

0 * T m,i( t ) < A < T c Vt, m e {l, . . . , m}, 1 e {l, . . . , l} (11) 

As a consequence of Eq. (11), the users appear synchronous at 
the receiver. 

Independent Ravleigh distributed multiplicative channel 
coefficients : We assume the channel coefficients a m/1 (t) in 
Eq. (8) between symbol intervals are modeled as independent 
zero-mean circular white Gaussian noise processes and remain 
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constant during each symbol period. In particular, for n,n' = 
...,«>, m,m' e {1,...,M} and 1, l'e { 1, . . . , L} . 

o^nTb+AJ-NfO, G^) (12) 

a m . 1 (nT b +A)= a m/1 (nT b +T c +A) = »• a m , x (nT b + (G-l ) T c +A) (13) 
5 E{a m , x (nT b +A)cv, r (n'T b +A) }=6 n , A,A,r^ A (14) 

where 5 k/k < denotes the Kronecker delta function, that is 5 k/k *=l, 

if k=k', and 5 k/k '=0 otherwise. 
Q Temporally white observation noise : The noise vector w(kT c +A) 

^ is modeled as a zero-mean circular white Gaussian noise 

lA) process with a known NxN covariance matrix R. For k = 
jJl . . . , °°, we have 

s" 

^ E{w(kT c +A) } = 0 Nxl 

;S E{w(kT c +A)w H (k'T c +A) } - R5 k/k - (15) 

FU where 0 Nxl denotes a Nxl vector of zeros. 

^5 Discrete-time state-space model : For notational convenience, 
we rewrite the kth sampled array output after beamforming, 
given by Eq. (10), as follows 

Y k = f k "r k + f k H W k (16) 

where y k £y ( kT c +A) , f k £f ( kT c +A) and w k £w(kT c +A). Here we use the 
20 subscript to denote the /cth sample. 

Let Xk, n ,i denote the attenuated received spectrum signal 
at time k for the mth user via the 1th path, i.e. 

Xk,»,i a b£ m c k/m a k/m/1 (17) 
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where ot krm<1 &a mrk (kT c +A) . Let 



IkJtLL Vk 6 {...,- 2(G - 1) ,-(G - 1) ,0, G - 1,2(G - 1),...} 
0. otherwise 



denote the realization of the attenuated and phase shifted 
spread spectrum Xk, m ,i every G samples. Since the product 

b k,m c k,m yields +1, and from our statistical assumptions on the 

channel fading, V k , m/1 is complex white Gaussian distributed 
random variable, such that 

2 

V k . m .! ~ N* 0 ' 0 ™,!*' for k 6 { - - -,-2(G-l) , - (G-l) , 0,G-1,2 (G-l) , . . . } (19) 



and U k , m ,i = 0, otherwise. Let 

[o, Vk 6 { . . . ,-2(G - 1) ,-(G - 1) ,0, G - 1,2(G - 1) , . . . } 

h k = I . . . . (20) 

II, otherwise 



The attenuated phase shifted received spread spectrum signal 

Xk, m ,i can be written in the following scalar recursive form 

Xk,m,l = hk C k,m C k-l,mXk-l,m, 1 + Vk,m, 1 (21) 

We group the received signals Xk,m,i at time index k in a MLxl 
vector Xk as follows 

H 

X k" ( Xk, 1, 1, . . . , Xk, 1,L • • • / Xk,M, 1, * • • t Xk,M,lJ ( 22 ) 
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, H 



If we define V k a (V k , 1#1 , . . . , V k , 1#L , . . . , V k/M , lf . . . , V k , M/L ) H , then we 



have the following vector state equation for the time 
evolution of the received spread spectrum signals for all user 
and all paths 



Z k * = A *Xk-i + v k 



(23) 



where 



v k - NfCW.Q) ,V k e { 
V k = 0, otherwise 



. ,-2(G-l) ,-(G-l) ,0,0-1,2(0-1) , 



.} (24; 
(25) 



and 



A k =diag 



h kPk-lA c W: * * ' h l«Fk-lA c k f l' * * * ' h IP k- 1^= k,M' • • • * h )^k- l^k.M 



LTirres 



L Tines 



(26) 



Q = diag (c^ v . . , a J^. . . , a^ v . . , a J fL ) (27) 

where. A k and Q are MLxML diagonal matrices. 

Let (J) denote the vector of all angle of arrivals for all 

M users. 

♦ H = (^lA' • • • ' <K,:l/ • • • ' *t* m,1' • • • * <t> M ,L) H ( 28 ) 

We use A^) to denote the N x ML matrix response of the array, 
defined as follows 
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A(t) = (a <*!,!) , . . . ,a(+ 1/L ) ,..wa(Vi)/.. • /a(t M , L ) H (2 9) 

The ( (m-1) L+l) the column vector in A(+) corresponds to the 
antenna array vector response of a signal impinging onto the 
antenna array from angle + m#1 . 
5 Using Eqs . (23) and (29), the observation Eq. (16) is 

equivalently written as 

Yk = f?A(4»X k + f£w k (30) 

Equations (23) and (30) from the discrete-time, state-space 
O version of the DS-CDMA beamforming linear antenna array 

ffD system. 

| Ja 

N 1 . Optimal Chip-Rate Switched-Beam Design 

y * 

|Jl In this subsection, we formulate the estimation 

^ objectives. The aim is to design the time varying switched- 

;D beam beamforming vector f k , such that the minimum estimation 

115 error in estimating the received spread spectrum signal Xk 
□ defined in Eq. (22) is obtained. We chose the beamforming 

vector f k to be a function of the data from time index 0 up to 
time index k-1, Furthermore, f k belongs to a set of fixed beam 
array patterns F. Thus, 
20 f k = f k (y 0 ,yi, • • .,y>c-iU F (31) 

for k = 1,2,.... The aim is to select the class beam patterns 
F and optimally select the sequence of beam-patterns from the 
set F, such that the estimation errors in estimating the 
spread spectrum signals are minimized. We choose the 
25 following optimization function 

h k = E{ ( Xk - Xk/k) H J(X* " Xwk) } (32) 
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where Xwk -E { Xk /Y k) is the conditional mean estimate of Xk/ 
conditioned on the data Y k ± <y 0 , y lt . . . , y k J J is some predetermined 
- user defined - weighting matrix. We choose the sequence of 

^ 

optimal beamf orming vector f 0 , f 1# f 2 , . . . ) , as follows 



q subject to Eq. (31). 

m l n the following subsection, we present several well 

u § 

% %J known beamf orming techniques. 

ill 2 . Beamf orming Techniques 

;10 We group our beamforming techniques according to 

1 SS3. 

• Fixed-beam system, where the antenna beams are fixed 
1U in azimuth. 

□ • Maximizing the response of the antenna array in the 

direction of a mobile unit. 
15 • Minimizing the average output power of the 

beamforming array, while ensuring unity response in 
a desired direction. 
We now present details on selecting the sets of beamforming 
array patterns based on the three beamforming techniques 
20 mentioned above. 

Fixed Beam System : The beamforming vectors in a fixed-beam 
antenna array system are designed to jointly span the entire 
area around the cell site. Let B denote the number of fixed 
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beams and ^ denote the ith angle of the beam where the antenna 
array is pointing. The set of fixed beams are given by 

F a = {at*!) , a(* 2 ) , . . a(t B ) } (34) 
Conventional Beam former : The beamforming vectors in a 
conventional beamformer are proportional to the steering 
vector in the known direction of the desired signal. The 
complex beamforming vector weights have equal magnitudes and 
the phases are set to the phases of the array vector response 
when the angle of arrival is in the desired direction 9. 

From the previous discussion, we assumed M x L steering 
vectors associated with M x L known angle-of -arrivals $ mtX . The 
set of conventional beamforming array vectors are given by 

F 2 = (0£(*i, i) • /0£(*i, l) • /0£(+ Mfl ) • /<x(4m,l) > < 35 ) 

Optimal Beamforming : The optimal beamforming vector is 
derived by minimizing the average output power of the 
beamforming array, while ensuring the unity response in the 
desired direction s 0 . Assuming f = f k for all k, we solve the 
following optimization problem 

H H 

min f E{yy } , subject to f s 0 = 1 (36) 

Using Lagrange's multiplier method, the optimal beamformer f 
can be expressed as follows 
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where D denotes the average signal, multi-user interference 
(MUI) plus noise power. From our statistical assumptions on 
our signal model, we can show the following result 

D = A(t)QA H (<h + R (38) 
5 Note that D is a positive definite matrix since Q and R are 
positive definite matrices. Thus, the set of optimal 
beamforming steering vectors are give by 

= [ p"^t($u> p'^tfrijj p~V<t>*a> 

° - t i i (39) 

a I \* Mf L)D" 3a (*M # L)J 

IT! 



111 III. OPTIMAL SWITCHED -BEAM DESIGN 

ffib We now summarize our state-space model, estimation 

objectives and present the optimal switched-beam sequence f k in 
the following theorem: 
Theorem 3.1 Gi ven 

System model : The multipath Rayleigh fading DS-CDMA spread 
15 spectrum communication system sampled at chip rate at the 
antenna array output after beamforming is described by the 
following discrete time equation 

X k+ i = A k+1 x k +V k+1 (40) 

y k = f?A<J) Xk + f k N w k ,forf k eF <4i> 
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where k is the time index, and Xk is the state vector 
representing the received attenuated and phase shifted 
versions of the DS-CDMA signals. The state transition matrix 
A k is given by Eq. (26) and the process noise v k is a zero-mean 
5 complex white Gaussian noise sequence give by Eqs. (24) and 

(25) . f k is the beamforming vector and it is an element of a 
finite set F of preselected steering vectors. w k is the noise 
vector at the antenna array, and w k is modeled as a zero-mean 
complex white Gaussian measurement noise sequence with 

O 

positive definite covariance R. Finally, we assume that the 
jfi process noise v k and observation noise w k are statistically 

!J1 uncorrelated, that is Efi^v^} = 0 for k,l = 0,1,2,... 

i"* Known Model joarameters : The angle-of-arrivals + for all spread 

w spectrum users and their associated steering vectors are 

1|5 known, thus A(+) is assumed known. The spreading sequences 

1 1 1. 

|«* also are known to the receiving antenna array, thus A k are 

known for all k. Finally, the noise covariance matrices Q and 
R are also known. When some of the model parameters are not 
known, then the unknown parameters can be estimated. Model 

20 parameter estimation is not investigated here. Standard 
techniques may be used for this purpose. 

Constraints : We chose F to be equal to one of the three 
beamforming vector sets discussed in Section II. 2. Only one 
beamforming vector f k eF may be used at any one time. The 
25 active beamforming vector at time k is selected via one of the 
following control laws. 
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F k = F(k, A((f>) , A k , Q, R) € F (42) 

f k = f(k, A((|)) , A k , Q, R, Y k-1 ) £ F (43) 

The closed-loop control function f ( k, A ( <h , A k , Q, R, Y k _ a ) - to be 
determined - at time k makes use of all available information 
up to time k, Y k denotes the measurement sequence 
(Yo/ Yi* • • / Yk) * Note that f k is a function of the data Y k _ x and 
thus the measurement y k is also a function of Y k _! • The open- 
loop control function f ( k, A ( ♦) , A k/ Q, R) at time k is a function 
of the known model parameters and a function of time. 
fObiective Functiqps : Compute the chip-rate beamf orming vector 

sequences F G-1 A (£ 0 , f 2 , . . . , f G _ 1 ) and F^Af^, f 1# . . . , f^) which 

respectively maximize the terminal cost function h^ = E { (Xg-i" 
Xg-iIg-i) H J(Xg-i-Xg-iIg-i) | as follows 

*g-i = ar^maxE {(z G _ x - Z G -i\g-i) H j (Zg-i " /g-i|g-i)> (44) 



F G _, = argmax^lC^.! - Z G -i\g-i) H J (Z G -i ~ Z G -x\g-x)) (4 5) 



Fa-x 



Then the optimal control sequences 
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( f 0 , • . . , f G -i) and ( f c • • • ' f G-i) obtained from (44) and (45), 

respectively, are computed off-line by solving a nonlinear 
deterministic optimization problem. In particular 

1 . f k = f k , for k = 0,1,... ,G - 1 



where P^b-! (F^UE { (Xg-i-Xg-iIg-i) Xg-i'Xg-iIg-i) "^g-i, F^) } denotes the 
error covariance matrix of Xg-i given the data Y^, P g -iIg-i ( f g-i) 

s e s 
: J c 

^ is the solution of the discrete-time Riccati equation, 

Y$ computed by running a Kalman filter, as follows 

h P k|k-l( F k-l> = Afc^-llk-l^k-l)*!! + Q (46) 

in 



S k!k-l< F k> = f?A((|»P ll|k _ 1 (F k _ 1 )A H (*)f k + f k Rf^ (47) 



for k = 0,...,G-1 and P-ii^F^) = 0^^, 

Proof : The proof follows similarly to the proof in L. Meier, 
15 J. Peschon and R.M. Dressier, Optimal Control of Measurement 
Subsystems, IEEE Transactions on Automatic Control, Vol. 12, 
No. 5, pp. 528-536, October 1967, where dynamic programming is 
used to compute the optimal control for controlling the 
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measurement equations in a Linear Quadratic Control (LQC) 
problem. 

Thus, the error covariance matrix P G -ib-i ( f g-i ) / as shown 
from the Riccati equations above, is data independent. 
Furthermore, the open-loop and closed-loop beam-forming vector 
sequences are identical. The optimal sequence is derived by 
using one of: 1) dynamic programming, 2 ) enumeration or 
3) enumeration via optimal pruning. 

Dynamic programming is not suitable for this particular 
application since the dynamic programming solution requires 
the quantization of a large dimensional error covariance 
matrix. The optimal beamforming sequence via enumeration is 
computed by considering all possible beamforming vector 
sequences and choosing the sequence that yields the minimum 
terminal cost. In A. Logothetis and A. Issakson, On Sensor 
Scheduling via Information Theoretic Criteria, American 
Control Conference, 1999, an algorithm that reduces the 
complexity in searching the optimal sequence is derived. The 
algorithm in that paper is an enumeration scheme with optimal 
pruning where sets of sequences are deleted without the 
possibility of removing the optimal sequence. 

IV. EXEMPLARY DEVICE 

As can be seen from the discussion above, the beam 
scheduling sequence does not require any of the observations 
to be determined. Further, it does not need to be revised at 
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a very high frequency since the environmental conditions that 
require revision of the beam scheduling sequence occur at a 
very slow rate compared to the processing speed of modern 
communications electronics. Accordingly, in one exemplary 
5 embodiment, the beamforming sequence can be computed offline 
and stored in a look-up table with the look up table being 
updated at regular intervals, such as once every frame. 

Figure 2 is a block diagram of the inventive switched 
beam antenna array receiver of the present invention. As 
fib shown, each of a plurality of antennas 112^ 112 2/ ..., 112 N 
|fj receives RF signals and outputs those signals on lines 113 lf 

113 2/ ..., 113 N to an analog N x 1 switched beamforming circuit 
■* s s 114. In accordance with the present invention as described in 

^ detail above, circuit 114 is an Nxl beamforming circuit that 

Jf^ generates a single beam at any given instant on line 115 in 

accordance with a schedule as controlled by a beamforming 
^ control circuit 140. 

The particular beam output on line 113 is switched at a 
very high rate, such as chip rate. As just noted, however, 
20 the switching schedule can remain unchanged for a relatively 

long period, e.g., several milliseconds. Accordingly, in this 
exemplary embodiment, the beam schedule is calculated and 
updated offline at, for instance, the frame rate of the 
communications system, and is stored in a look up table in a 
25 memory 142. The control circuit 142 retrieves the beamforming 
schedule from the look up table to control circuit 114. 
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In accordance with the present invention, the beamf orming 
control circuit also needs direction of arrival estimates of 
the various signals in the reception region (cell) to properly 
set the beam schedule. Accordingly, a direction of arrival 
5 estimator circuit 144 estimates the direction of arrival of 
the various user signals and provides those estimates to the 
beamforming control circuit 140. Direction of arrival 
estimators are well known in the prior art and therefore are 
not described in detail herein. 

ifep The single beam output from circuit 144 onto line 115 is 

in 
■ =• 

frequency down converted by frequency downconverting circuit 

^ 116 and converted from analog-to-digital by an analog-to- 

^} digital converter 118. The digital beam signal is input to a 

;^ minimum variance flat-fading multipath and multibeam channel 

J5 estimator 122 to generate L path estimates and L path estimate 

fy errors for each of M simultaneous transmitters. Suitable 

.'tS3. 

minimum variance flat-fading multipath and multi-user channel 
estimator circuits and algorithms are well known in the prior 
art. The invention is not limited to use in flat fading 
20 channels. The inventive algorithm may be implemented by any 
form of digital processing apparatus. For exemplary purposes 
in connection with Figure 2, the algorithm may be considered 
to be performed by a DSP in block 122. 
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V. CONCLUSION 

In this specification, we considered the problem of 
optimally scheduling a multibeam antenna array for DS-CDMA 
communication systems operating in a multipath Rayleigh flat- 
fading channel. We derived a scheme that yields the minimum 
estimates of the received spread spectrum signals by computing 
optimal chip-rate beamforming scheduling (switching) 
sequences. More particularly, we formulated the received 
communication system as a linear, Gaussian, time-varying 
state-space model, the individual paths of the individual 
transmitted signals in the channel defining the "states" of 
the state-space model. We performed minimum mean square error 
estimation of the state of the state space model (i.e., the 
mean square estimation of individual paths of the individual 
signals in the received spread spectrum channel) to generate a 
beamforming schedule with optimal signal/interference-plus- 
noise ratio. This is a traditional use of Kalman filtering, 
namely, deriving the individual states from a signal 
containing a linear combination of states. 

The advantages of our proposed scheme are: 1) only one 
beam is generated at each time instant (the generated antenna 
beam is used to estimate all spread spectrum users and their 
transmission paths) ; 2) the beam pattern of the antenna array 
is selected among a predetermined set of beam patterns; 3) the 
beam pattern changes at chip rate; 4) the switching sequence 
of the antenna beams is data independent when the direction of 
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arrival (DOA) of the received spread spectrum signals are 
known (thus, the switching sequence is determined a priori) ; 
and 5) the optimal switched beam sequence is designed to yield 
minimum estimates of the communication channel. 

Having thus described a few particular embodiments of the 
invention, various alterations, modifications, and 
improvements will readily occur to those skilled in the art. 
Such alterations, modifications and improvements as are made 
obvious by this disclosure are intended to be part of this 
description though not expressly stated herein, and are 
intended to be within the spirit and scope of the invention. 
Accordingly, the foregoing description is by way of example 
only, and not limiting. The invention is limited only as 
defined in the following claims and equivalents thereto. 



